Dexamethasone versus methylprednisolone for multiple organ dysfunction in COVID-19 critically ill patients: a multicenter propensity score matching study

Background Dexamethasone usually recommended for patients with severe coronavirus disease 2019 (COVID-19) to reduce short-term mortality. However, it is uncertain if another corticosteroid, such as methylprednisolone, may be utilized to obtain better clinical outcome. This study assessed dexamethasone’s clinical and safety outcomes compared to methylprednisolone. Methods A multicenter, retrospective cohort study was conducted between March 01, 2020, and July 31, 2021. It included adult COVID-19 patients who were initiated on either dexamethasone or methylprednisolone therapy within 24 h of intensive care unit (ICU) admission. The primary outcome was the progression of multiple organ dysfunction score (MODS) on day three of ICU admission. Propensity score (PS) matching was used (1:3 ratio) based on the patient’s age and MODS within 24 h of ICU admission. Results After Propensity Score (PS) matching, 264 patients were included; 198 received dexamethasone, while 66 patients received methylprednisolone within 24 h of ICU admission. In regression analysis, patients who received methylprednisolone had a higher MODS on day three of ICU admission than those who received dexamethasone (beta coefficient: 0.17 (95% CI 0.02, 0.32), P = 0.03). Moreover, hospital-acquired infection was higher in the methylprednisolone group (OR 2.17, 95% CI 1.01, 4.66; p = 0.04). On the other hand, the 30-day and the in-hospital mortality were not statistically significant different between the two groups. Conclusion Dexamethasone showed a lower MODS on day three of ICU admission compared to methylprednisolone, with no statistically significant difference in mortality. Supplementary Information The online version contains supplementary material available at 10.1186/s12879-024-09056-y.

Corticosteroids have anti-inflammatory and immunosuppressive actions that reduce the production of pro-inflammatory cytokines such as IL-6 [15].Previous studies have shown improved prognosis and survival rate in patients with moderate to severe COVID-19 using corticosteroids therapy, including dexamethasone, methylprednisolone, and hydrocortisone [10,14,[16][17][18].However, the use of hydrocortisone or methylprednisolone in this patients population was not supported by solid evidence [19].The RECOVERY trial is a randomized control study that included 6425 hospitalized patients with COVID-19.The study found that patients using dexamethasone along with standard therapy had a lower mortality rate at 22.9% compared to patients who received standard treatment alone at 25.7% (rate ratio, 0.83; 95% confidence interval (CI), 0.75 to 0.93; P < 0.001) [19].
Methylprednisolone had greater lung tissue-to plasma penetration in rats when compared to dexamethasone.This effect of methylprednisolone might be beneficial to reduce lung damage and prevent respiratory complications [20].Studies investigating methylprednisolone in patients with severe COVID-19 showed improved survival [10].Another randomized controlled trial (RCT), including 68 hospitalized patients admitted to the ICU with COVID-19, compared pulse methylprednisolone to the standard of care.This study found the administration of pulse methylprednisolone at the early phase of illness significantly increased survival compared to standard of care (HR 0.293; 95% CI 0.154-0.556;p < 0.001).The use of methylprednisolone helped improve pulmonary symptoms and inflammatory markers [16].
Both intravenous methylprednisolone and dexamethasone, administered in a 2 mg/kg/day and 6 mg/day dose, respectively, demonstrated similar efficacy in reducing hospital length of stay and oxygen therapy needs in patients with severe viral infection manifestation [21].In contrast, a small RCT that enrolled 86 critically ill ventilated patients with COVID-19; found that patients receiving methylprednisolone have shorter hospitalization and less dependency on mechanical ventilation (MV) compared to dexamethasone [22].In non-critically ill COVID-19 patients, pneumonia treated with high-dose methylprednisolone showed lower C-reactive protein (CRP), and D-dimer compared to high-dose dexamethasone [1].
Multiple organ dysfunction score (MODS) is a validated tool that measures the severity of MOD in the intensive care unit.MODS consisted of 24 scores that could be used as a prognostic indicator if calculated at ICU admission or as an outcome once calculated during ICU stay.Higher MODS are indicative of a higher probability of ICU mortality [23].During the pandemic, severe COVID-19 infection was associated with MOD involving cardiovascular, renal, and pulmonary systems [24].The involvement of many vital organs in patients with coronavirus infection made the MODS as the most relevant indicator for a patient's complex condition.National Institutes of Health Guidelines recommended using systemic corticosteroids, including dexamethasone and methylprednisolone for hospitalized patients with COVID-19 on supplementary oxygen [19].Most previous studies investigated the survival benefits and the need for MV.Corticosteroids reduced the risk of endotracheal intubation as indicative of improved lung function [25].However, the corticosteroid of choice for better clinical outcomes in critically ill patients with COVID-19 remains undetermined.Thus, we conducted this study to assess the short termclinical and safety outcomes of dexamethasone compared to methylprednisolone in critically ill patients with COVID-19.

Study design
This research is a component of the Saudi Critical Care Pharmacy Research (SCAPE) platform, which has conducted numerous investigations into the safety and efficacy of diverse treatments and therapies for critically ill patients (Saudi Critical Care Pharmacy Research (SCAPE), 2023).The present study's design and methodologies draw upon those employed in our previously published research [6].This is a multicenter retrospective cohort study.It included adult critically ill patients with confirmed COVID-19 admitted to the ICUs between March 01, 2020, and July 31, 2021.COVID-19 was diagnosed using Reverse Transcriptase-Polymerase Chain Reaction (RT-PCR) nasopharyngeal or throat swabs.We categorized the included patients into two groups depending on the type of corticosteroids used within 24 h of ICU admission (Dexamethasone or Methylprednisolone).We followed eligible patients until hospital discharge or death.The King Abdullah International Medical Research Center (KAIMRC) authorized the study in February 2022 (Ref.#NRC22R-074-02).Due to the study's retrospective observational nature, informed consent from study participants was waived.All methods were performed following relevant guidelines and regulations.

Study participants
All adult patients (age ≥ 18 years) with confirmed COVID-19 admitted to the ICUs during the study period who received either dexamethasone or methylprednisolone were assessed for eligibility.Patients were excluded if they received dexamethasone or methylprednisolone before ICU admission or those who were initiated on corticosteroids more than 24 h after ICU admission.Other exclusion criteria include patients who used corticosteroids as concomitant or sequential therapy, ICU discharge, or death within the first 24 h of ICU admission, and those labeled as "Do-Not-Resuscitate" (Fig. 1).

Study setting
This study was conducted at five hospitals in Saudi Arabia: King Abdulaziz Medical City (Riyadh and Jeddah), King Abdulaziz University Hospital (Jeddah), King Abdullah bin Abdulaziz University Hospital (Riyadh), and King Salman Specialist Hospital (Hail).The selection of these centers was based on the geographic distribution, availability of electronic records, and the center's willingness to participate in the national project.The primary site for this multicenter study was King Abdulaziz Medical City (KAMC-Riyadh), a tertiary care center.

Data collection
Each patients' data were collected and managed using Research Electronic Data Capture (REDCap ® ) software hosted by KAIMRC.We gathered patients' demographic data, comorbidities, vital signs, laboratory tests, baseline severity scores (i.e., MODS, Acute Physiology and Chronic Health Evaluation II (APACHE II), and Sequential Organ Failure Assessment (SOFA)), renal profile, acute kidney injury, use of prone positioning, MV use and parameters (e.g., PaO2/FiO2 ratio, FiO2 requirement) within 24 h of ICU admission.The Glasgow Coma Score (GCS) was collected in the study, with sedation assessment using the Richmond Agitation-Sedation Scale (RASS) for mechanically ventilated (MV) patients.Conscious level evaluation for non-MV patients or those on non-invasive MV utilized the GCS score.Verbal scores were excluded in invasively intubated and sedated patients due to their inability to vocalize, and assessment focused on eye opening and motor scores, allowing for the computation of the GCS sum score in intubated patients.For noninvasive MV patients, the best GCS was used for conscious level assessment, and for invasive MV patients, GCS scores were measured during sedation vacations, offering a comprehensive understanding of their neurological status.We gathered information on liver function tests (LFTs), coagulation profile (i.e., INR, aPTT, fibrinogen, D-dimer), and inflammatory markers (Ferritin, procalcitonin, and creatine phosphokinase (CPK)) within 24 h of ICU admission.Data on the use of corticosteroids therapy (type, dose, and duration), the timing of corticosteroids initiation, and tocilizumab use were collected for the eligible patients.

Outcomes
The primary outcome was the multiple organ dysfunction on day three of ICU admission.The secondary outcomes were mortality, hospital Length of Stay, ICU Length of Stay, respiratory failure requiring MV, ventilator-free days (VFDs) at 30 days, and ICU-acquired complications (new-onset atrial fibrillation, AKI, liver injury, hospital-acquired infection, pneumonia, and secondary fungal infection).Outcomes definitions are shown in the Additional file 1.

Statistical analysis
We presented continuous vaiables as mean and standard deviation (SD), or median and lower quartile (Q1) and upper quartile (Q3), as appropriate and categorical variables as number (percentage).The normality assumptions were assessed for all numerical variables using a statistical test (i.e., Shapiro-Wilk test) and graphical representation (i.e., histograms and Q-Q plots).Model fit was assessed using the Hosmer-Lemeshow goodness-of-fit test.
Baseline characteristics and outcome variables were compared between the two study groups.We used the Chi-square or Fisher's exact test for categorical variables.We compared the normally distributed continuous variables using the student t-test and other non-normally distributed continuous variables with the Mann-Whitney U test.Multivariable Cox proportional hazards regression analyses were performed for the 30-day and in-hospital mortality.The proportionality assumption was assessed before fitting the cox model.Visual assessment was performed to evaluate the hypothesis by plotting a log(-log) plot and testing the correlation of scaled Schoenfeld residuals with rank-ordered time.
Multivariable logistic regression analysis and negative binomial regression were used for the other outcomes considered in this study.Regression analysis was done for the study outcomes after considering propensity scores (PS) as covariates in the model.These scores are generated through PS analysis after considering all relevant covariates, which includes patient's age and MODS within 24 h of ICU admission.The odds ratios (OR), hazard ratio (HR), or estimates with the 95% CI were reported as appropriate.No imputation was made for missing data as the cohort of patients in our study was not derived from random selection.
The Proc PS match procedure (SAS Cary, NC) was used to match patients who received methylprednisolone therapy (active group) to patients who received dexamethasone therapy (control group) (1:3 ratio) based on patient's age and MODS within 24 h of ICU admission.A greedy nearest neighbor matching method was utilized; one patient who received methylprednisolone therapy (active) paired with three patients who received dexamethasone (control group), which eventually produced the smallest within-pair difference among all available pairs with treated patients.The difference in the logits of the PS for pairs of patients from the two groups was matched if it was less than or equal to 0.5 times the pooled estimate of the standard deviation.We considered a P value of < 0.05 statistically significant and used SAS version 9.4 for all statistical analyses.

Results
Among 1385 patients who were screened during the study period, 526 patients enrolled after applying the study exclusion criteria.Of these patients, 455 were given dexamethasone, while 71 were given methylprednisolone therapy 24 h after ICU admission.After PS matching (1:3), 264 patients were included (66 patients received methylprednisolone, while 198 patients received dexamethasone) according to the selected criteria.The median dose of methylprednisolone and dexamethasone per day was 80 mg (80, 120) and 6 mg (6, 6), respectively.Methylprednisolone was administered for a median duration of six days (4.7, 10) compared with ten days (9,11) for the dexamethasone group.

Demographic and clinical characteristics
Most of the patients were male (65.8%) in the whole cohort, and the average age was 62.0 ± 15.25 years.Diabetes mellitus (60.0%) was the most common comorbidity, followed by hypertension (56.3%), dyslipidemia (21.9%), and chronic kidney disease (11.2%).There was a statistically significant difference in the baseline characteristics between the two groups before PS matching.
Patients who received dexamethasone have a lower PaO2/FiO2 ratio, higher use of inotropes/vasopressors, and GCS at admission than methylprednisolone.After PS matching based on the patient's age and MODS within 24 h of ICU admission, most of the baseline demographics and clinical characteristics were well balanced between the two groups, except for the higher total bilirubin, INR, and P/F ratio in the methylprednisolone group.There were no significant differences in severity scores within 24 h of admission (APACHE II, SOFA scores, and MODS) between the two groups, regardless of PS matching.In addition, the proportion of tocilizumab use within 24 h of ICU admission and the concomitant use of nephrotoxic medications during ICU stay were similar between the two groups (Table 1).

Multiple organ dysfunction (MODS)
Details of the MODS on day three are shown in (Table 2).Best GCS, CVP, heart rate, pressure adjusted heart rate, and PaO2/FiO2 were significantly different between the two groups.The PaO2/FiO2 on day three of ICU admission was higher in the methylprednisolone group than dexamethasone; however, it was higher at baseline (Table 1).In crude analysis, the MODS on day three of ICU was higher in patients who received methylprednisolone when compared to dexamethasone.Similarly, at regression analysis, patients in the methylprednisolone group have a higher MOD score (beta coefficient: 0.17 (95% CI 0.02, 0.32), P = 0.03) (Table 3).

Complications during ICU stay
Among the non-MV patients within 24 h of admission, patients who received methylprednisolone have a higher odd of respiratory failure that requires MV; however, it was not statistically significant (OR 1.60, 95% CI 0.53, 4.85; p = 0.41).Patients who received methylprednisolone had a higher odds of infection complications during ICU than dexamethasone; however, didn't reach a statistically significant difference.Moreover, other complications during the stay, such as new-onset atrial fibrillation, AKI, and liver injury, were similar between the two groups (Fig. 2a, b).

30-day and in-hospital mortality
The prevalence of 30-day mortality in the crude analysis was 34.8% in the dexamethasone group compared with 41.7% in patients who received methylprednisolone, which was not statistically significant (P = 0.34).Similarly, multivariable time dependent cox proportional hazards regression analysis, revealed a similar 30-day and inhospital mortality (HR 0.997, 95% CI 0.97, 1.03; p = 0.85 and HR 1.02 95% CI 1.0, 1.04; p = 0.053 respectively) (Table 3).

Discussion
This retrospective multicenter study compared clinical and safety outcomes between dexamethasone and methylprednisolone among 264 critically ill patients with COVID-19.After PS matching and linear regression model analysis, we found a statistically significant lower MODS among patients treated with dexamethasone compared to methylprednisolone.Similar findings were reported after adjusting for age and MODS within 24 h of ICU admission.Indeed, our study found that patients treated with methylprednisolone had significantly lower GCS, higher PaO 2 /FiO 2 , and pressure-adjusted heart rate.One possible explanation for higher PaO 2 /FiO 2 is the well-established higher lung tissue-to-plasma ratio of methylprednisolone [26][27][28].Although is it controversial, it has been considered that dexamethasone penetrates the blood-brain barrier relatively more than other corticosteroids [29][30][31].This could explain the significantly higher GCS scores among patients treated with dexamethasone in our study.It is important to note that GCS assessment in invasive MV is limited by the uncertainty surrounding GCS assessments in the context of intubation.The CoDEX trial investigated the effect of dexamethasone on COVID-19 patients and concluded that dexamethasone resulted in a significantly higher number of free-from MV days than the standard of care alone despite the similar PaO 2 /FiO 2 ratios of both groups [32].Furthermore, in the RECOVERY trial, dexamethasonetreated patients had a lower need for renal-replacement therapy and invasive MV [16].Despite using a lower dexamethasone dose compared to CoDEX trial and having a higher CRP and lower PaO2/FiO2 ratio at admission, our dexamethasone-treated patients had lower MODS [32].
Regarding secondary outcomes in our study, our study showed a higher odds of infection complications during ICU among methylprednisolone-treated cohort than dexamethasone; however, didn't reach to a statically significant difference.A large meta-analysis suggests higher secondary infections following corticosteroid use among influenza patients [33].However, similar findings were not observed in other studies that mainly looked at COVID-19 patients [34].The Metcovid trial included suspected and confirmed cases of COVID-19 and investigated the effect of methylprednisolone on mortality [35].The similar rates of hospital-acquired infection in the Metcovid study could be explained by the absence of use of other immunomodulator COVID-19 therapies.In our study, a major factor could have contributed to the higher hospital-acquired infection rate with     Patients who received either Enoxaparin 40 mg daily or UFH 5000 Unit three times daily were grouped under the "standard dose VTE prophylaxis.Any patient who received higher than standard dose but not as treatment dose (Enoxaparin 1 mg/kg q12hr or 1.5 mg/kg q24hr or UFH infusion) was categorized as receiving "High VTE prophylaxis dose".On the other hand, lower VTE prophylaxis considered for patient who received Enoxaparin < 40 mg/day or Unfractionated heparin (UFH) < 5000 Units three times daily/day) methylprednisolone.Methylprednisolone is a relatively more potent immunosuppressant than dexamethasone [35].Although it was statistically non-significant, patients treated with methylprednisolone had higher ICU and hospital lengths of stay.Several studies showed longer hospital and ICU lengths of stay result in a higher infection rate [36,37].Most importantly, the higher MODS and infection rates were not translated into a higher mortality rate with methylprednisolone than dexamethasone.While we did not statistically adjust for SOFA scores or age, the relatively young patient population and low predicted mortality rate based on SOFA scores among our patients could have impacted the mortality rates in our analysis.Similarly, other studies demonstrated a reduction in 28-days and in-hospital mortality [10,[12][13][14].Three of the studies mentioned included more elderly patients and had used more immunomodulators [10,12,13].Another study included patients with much lower CRP than ours, indicating a less acute COVID-19 inflammation [14].These factors could be attributed to the lower mortality rates in other studies.In contrast to other studies, we reported similar outcomes between methylprednisolone and dexamethasone on ventilator-free days and respiratory failure requiring MV [10,12,34].The low need for MV at admission and lower ICU admission rates in the previously published studies might explain the differences in our findings.

Table 1 Summary of demography and baseline characteristics
Our study has a few limitations, including its retrospective nature which limit capturing some important variables such a; medications used prior to ICU admission, as well as onset of complications such as infections and relatively small sample size.In addition, our primary outcome MODS was only assessed at day three which limits the evaluation of the variation in MODS scores between the two groups over a continuous period.On the other hand, this study has some strengths, including multicenter patient recruitment, and PS matching was used to eliminate a greater portion of bias and create a balanced dataset.

Conclusion
This study showed a lower MODS on day three of ICU admission in critically ill patients with COVID-19 treated with dexamethasone.However, both interventions had a similar mortality rate.Future RCTs are needed to confirm our findings and to assess MODS over an extended period during ICU stay.

Fig. 1
Fig. 1 Eligibility criteria flowchart sum test is used to calculate the P-value ^^ Chi square ** Fisher's exact teat is used to calculate P-value a Nephrotoxic medications/ material included IV Vancomycin, Gentamicin, Amikacin, Contrast, Colistin, Furosemide, and/or Sulfamethoxazole/trimethoprim b

Fig. 2 a
Fig. 2 a Complications during ICU stay.b Infection complications during ICUs

Table 2
Details of the multi-organ dysfunction on day 3 Wilcoxon rank sum test is used to calculate the P-value * T -Test ^

Table 3
The outcomes of critically ill patients with COVID-19 after propensity score matching ^ Wilcoxon rank sum test is used to calculate the P-value ^^ Chi-square test is used to calculate the P-value $ Time dependent cox proportional hazards regression analysis used to calculate HR and p-value $* Generalized linear model is used to calculate estimates and p-value